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Abstract 

We invite the reader (presumably an upper level undergraduate 
student) to a journey leading from the continent of Classical Mechanics 
to the new territories of Quantum Mechanics. We'll be riding the 
symplectic camel and have William of Occam as travel companion, so 
no excess baggage is allowed. The first part of our trip takes us from the 
symplectic egg to Gromov's non-squeezing theorem and its dynamical 
interpretation. The second part leads us to a symplectic formulation 
of the quantum uncertainty principle, which opens the way to new 
discoveries. 



Prologue 

What's in a name? That which we call a rose by any other 
name would smell as sweet. 

Romeo and Juliet, Act 2, Scene 2 (W. Shakespeare) 

Take an egg -preferably a hard boiled one, and cut it in half along its 
middle using a very sharp knife. The surface of section will be roughly 
circular and have area irr 2 . Next, take a new egg of same size, and cut it 
this time along a line joining the egg's tops, again as shown in Figl. This 
time we get an elliptic surface of section with area irR 2 larger than that 
of the disk we got previously. So far, so good. But if you now take two 
symplectic eggs, and do the same thing, then both sections will have exactly 
same area! Even "worse", no matter along which plane passing through 
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the center of the egg you cut, you will always get sections having the same 
area! This is admittedly a very strange property, which you probably never 
have experienced (at least in a direct way) in everyday life. But what is a 
symplectic egg? The eggs we are cutting are metaphors for ellipsoids; an 
ellipsoid is a round ball that has been deformed by a linear transformation 
of space, i.e. a transformation preserving the alignment of three, or more, 
points. In mathematics such transformations are represented by matrices. 
Thus the datum of an ellipsoid is the same thing as the datum of a ball and 
of a matrix. What we call a symplectic egg is an ellipsoid corresponding 
to the case where the matrix is symplectic (we'll define the concept in a 
moment). The reason for which the only symplectic egg you have seen on 
your breakfast table is flat -a fried egg!- is that the number of rows and 
columns of a symplectic matrix must always be even. Since we are unable 
to visualize things in dimension three or more, the only symplectic eggs 
that are accessible to our perception are two dimensional. But what is a 
symplectic matrix? In the case of smallest dimension two, a matrix 



In higher dimensions, 4, 6, 8, etc. there are many more conditions: for 
instance 10 if the dimension is 4, 21 if it is 6, and n(2n + 1) if it is In. We 
will write these conditions explicitly in section 12. 1L 

So far, so good. But where do symplectic eggs come from, and what are 
they good for? Let me first tell you where symplectic matrices come from. 
They initially come from the study of motion of celestial bodies, which is 
really rich in mathematical concepts, some of these going back to the obser- 
vations of Tycho Brahe, and the work of Galileo Galilei and Johannes Kepler 
(these were the "Giants" on the shoulder's of which Isaac Newton stood). 
But the notion of symplectic matrix, or more generally that of symplectic 
transformation, did really have a long time to wait until it appeared explic- 
itly and was recognized as a fundamental concept. It was implicit in the 
work of Hamilton and Lagrange on classical and celestial mechanics, until 
the word "symplectic" was finally coined by the mathematician Hermann 
Weyl in his book The Classical Groups, edited in 1939, just before World 
War II. But still then, as Ian Stewart [30] reminds us, it was a rather baffling 
oddity which presumably existed for some purpose -but which? It was only 





is symplectic if it has determinant one: 



ad — be = 1. 



(2) 
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later agreed that the purpose of symplectic transformations is dynamics, 
that is the study of motion. Let me explain this a little bit more in detail: 
if we have a physical system consisting of "particles" (sand corns, planets, 
spacecraft, or quarks) it is economical from both a notational and compu- 
tational point of view to describe their motion (that is, their instantaneous 
location and velocity) by specifying a phase space vector, which is a matrix 
having only one column. For instance, if we are dealing with one single 
particle with coordinates (x,y,z) and momentum (p x ,Py,Pz) (the momen- 
tum of a particle is just its velocity multiplied by its mass m) the phase 
space vector will be the column vector whose entries are (x,y, z,p x ,p y ,p z ) 
If we have a large number N of particles with coordinates (xi,yi,Zi) and 
momenta (Px^Py^PzJ the phase space vector will be obtained by first writ- 
ing all the position coordinates and thereafter the momentum coordinates 
in corresponding order, their momenta. These vectors form the phase space 
of our system of particles. It turns out that the knowledge of a certain 
function, the Hamiltonian (or energy) function, allows us to both predict 
and retrodict the motion of our particles; this is done by solving (exactly, or 
numerically) the Hamilton equations of motion, which are in the case n = 1 
given by 

dx dH dp dH 

~dt = ~dp~ ' dt = ~~dx~' ( ' 

Mathematically these equations are just a fancy way to write Newton's sec- 
ond law F = ma. That is, knowing exactly the positions and the momenta 
at some initial time, we are able to know what these are going to be at any 
future time (we can actually also calculate what they were in the past). The 
surprising, and for us very welcome fact is that the transformation which 
takes the initial configuration to the final configuration is always a symplec- 
tic transformation! These act on the phase vectors, and once this action 
is known, we can determine the future of the whole system of particles, 
and this at any moment (mathematicians would say we are in presence of a 
"phase space flow"). The relation between symplectic transformations and 
symplectic matrices is that we can associate a symplectic matrix to every 
symplectic transformation: it is just the Jacobian matrix of that transfor- 
mation. In the simplest cases, for instance when no external forces act on 
the particles, these matrices are themselves the symplectic transformations. 

The symplectic egg is a special case a deep mathematical theorem dis- 
covered in 1985 by the mathematician Gromov [17), who won the Abel Prize 
in 2010 for his discovery (the Abel Prize is the "true" substitute for the No- 
bel Prize in mathematics, as opposed to the Fields medal, which is intended 
to mathematicians under 40). Gromov's theorem is nicknamed the "prin- 
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ciple of the symplectic camel" [31 [30], and it tells us that it impossible 
to squeeze a symplectic egg through a hole in a plane of "conjugate coordi- 
nates" if its radius is larger than that of the hole. That one can do that with 
an ordinary (uncooked) egg is easy to demonstrate in your kitchen: put it 
into a cup of vinegar (Coca Cola will do as well) during 24 hours. You will 
then be able to squeeze that egg through the neck of a bottle without any 
effort! 

The marvelous thing with the symplectic egg is that it contains quantum 
mechanics in a nutshell... er ... an eggshell! Choose for radius the square 
root of Planck's constant h divided by 2tt. Then each surface of section will 
have radius of h/2. In [BJ EJ El [14] I have called such a tiny symplectic egg 
a quantum blob. It is possible -and in fact quite easy if you know the rules 
of the game- to show that this is equivalent to the uncertainty principle of 
quantum mechanics. The thing to remember here is that a classical property 
(i.e. a property involving usual motions, as that of planets for instance), 
here symbolized by the symplectic egg, contains as an imprint quantum 
mechanics! The analogy between "classical" and "quantum" can actually 
be pushed much further, as I have shown with Basil Hiley [15]. But this, 
together with the notion of emergence |11| . is another story. 

Some of the ideas presented here are found in our Physics Reports paper 
|16j with F. Luef; they are developed and completed here in a different way 
more accessible to a general audience. 

1 Notation and terminology 

Position and moment vectors will be written as column vectors 



where the superscript indicates transposition. The integer n is unspeci- 
fied; we will call it the number of degrees of freedom. If the vector (x,p) T 
denotes the phase vector of a system of TV" particles, then n = 3N and the 
numbers xi,X2,X3, (resp. Pi,P2,P3) can be identified with the positions 




and the corresponding phase vector is thus 



Q = (^) T = c 



X\ , . . . , X n , pi , . . . , Pn) 
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x,y,z (resp. the momenta p x ,Py,Pz) of the first particle, £4,3:5, xq, (resp. 
P4,P5,Pe) with those of the second particle, and so on. This is not the only 
possible convention, but our choice has the advantage of making formulas in- 
volving symplectic matrices particularly simple and tractable. For instance, 

the "standard symplectic matrix" is here J = ( f 1 ) where 1^ is the 



n x n identity matrix and the n x n zero matrix. Note that 



J 2 = -I d , J T = J- 1 = -J. (4) 

2 The Symplectic Egg 

2.1 Symplectic matrices 

Let S be a (real) matrix of size 2n. We say that S is a symplectic matrix if 
it satisfies the condition 

S T JS = J (5) 

Clearly the standard symplectic matrix J is itself a symplectic matrix. 
Assume that we write the matrix S in block form 

where A, B, C, D are matrices of size n. It is a simple exercise in matrix alge- 
bra to show that condition (j^J) is equivalent to the the following constraints 
on the blocks A, B, C, D 

A T C = C T A , B T D = D T B , and A T D - C T B = I d . (7) 

Notice that the two first conditions mean that both A T C and B T D are 
symmetric). Observe that these conditions collapse to (j2J) when n = 1: 
in this case A, B, C, D are the numbers a, b, c, d so that A T C = ac and 
B T D = bd are automatically symmetric; the condition A T D — C T B = 
reduces to ad — be = 1 . 

The product of two symplectic matrices is a symplectic matrix: if S 
and S' satisfy © then (SS' fJSS' = S' T (S T JS)S' = S' T JS' = J. Also, 
symplectic matrices are invertible, and their inverses are symplectic as well: 
first, take the determinant of both sides of S T JS = J we get det(S T JS) = 
det J; since det J = 1 this is (detS") 2 = 1 hence S is indeed invertible. 
Knowing this, we rewrite S T JS = J as JS = (S~ 1 ) T J, from which follows 
that (S'~ 1 ) T JS'~ 1 = JSS^ 1 = J hence S 1-1 is symplectic. The symplectic 
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matrices of same size thus form a group, called the symplectic group and 
denoted by Sp(2n). An interesting property is that the symplectic group is 
closed under transposition: if 5 is a symplectic matrix, then so is S T (to 
see this, just take the inverse of the equality (S~ l ) T JS -1 = J). Since this 
means that a matrix is symplectic if and only if its transpose is, inserting 
S T in © and noting that (S T ) T = S we get the condition 

SJS T = J. (8) 

(A B\ ( A T C T \ 

Replacing S = ( ^ ^ J with S T = ( ^ T J the conditions ([7]) are thus 

equivalent to the set of conditions: 

AB T = BA T , CD T = DC T , AD T - BC T = I d . (9) 

One can obtain other equivalent sets of conditions by using the fact that 
5 _1 and (S~ r ) T are symplectic (see [8]). 

It is very interesting to note that the inverse of a symplectic matrix is 

^=(_ D c T r "if). (10) 

It is interesting because this formula is very similar to that giving the inverse 

] of a 2 x 2 matrix \ \\ with determinant one. The inversion 
-c a J \c a J 

formula fllOf) suggests that in a sense symplectic matrices try very hard to 

mimic the behavior of 2 x 2 matrices. We will see that this is actually 

the essence of symplectic geometry, and at the origin of the symplectic egg 

property! 

A last property of symplectic matrices: recall that when we wanted to 
show that a symplectic matrix always is invertible, we established the iden- 
tity (det S) 2 = 1. From this follows that the determinant of a symplectic 
matrix is a priori either 1 or —1. It turns out -but there is no really elemen- 
tary proof of this- that we always have detS" = 1 (see for instance §2.1.1 
in [8] where I give one proof of this property; Mackey and Mackey's online 
paper [20] give a nice discussion of several distinct methods for proving that 
symplectic matrices have determinant one. 

Conversely, it is not true that any 2n x In matrix with determinant one 
is symplectic when n > 1. Consider for instance 
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where o / 0; this matrix trivially has determinant one, but the condition 
AD T — BC T = 1^ in ([9]) is clearly violated unless a = ±1. Another simple 
example is provided by 



M 



R(a) 
R(P) 



where R(a) and R((3) are rotation matrices with angles a ^ f3 (this coun- 
terexample generalizes to an arbitrary number 2n of phase space dimen- 
sions) . 

2.2 The first Poincare invariant 

In what follows j(t), < t < 2ir, is a loop in phase space: we have ^(t) = 
x(t)\ 

^1 where x(0) = x(2tt), p(0) = p{2tt); the functions x(t) and p(t) are 

supposed to be continuously differentiable. By definition, the first Poincare 
invariant associated to j(t) is the integral 

J( 7 ) = & p dx= / p{tfx{t)dt. (12) 

J-y Jo 

The fundamental property -from which almost everything else in this paper 
stems- is that 1(7) is a symplectic invariant. By this we mean that if we 
replace the loop 7(f) by the a new loop S~f(t) where S is a symplectic matrix, 
the first Poincare invariant will keep the same value: I(Sj) = 1(7), that is 



pdx = d> pdx. (13) 
7 J Si 

The proof is not very difficult if we carefully use the relations characterizing 
symplectic matrices (see Arnol'd |2J, §44, p. 239, for a shorter but more 
abstract proof). We will first need a differentiation rule for vector- valued 
functions, generalizing the product formula d{uv)/dt = u(dv/dt) +v(du/dt) 
from elementary calculus. Suppose that 



u(t) 



are vectors depending on the variable t and such that each component Uj(t), 
Vj(t) is differentiable. Let M be a symmetric matrix of size n and consider 



/ U1 (t)\ 


, v(t) = 


Mt)\ 






Wit) J 
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the real- valued function u(t) T Mv(t). That function is differentiable as well 
and its derivative is given by the formula 



d_ 
~dt 



[u(t) T Mv{t)} = ii(t) T Mv(t) + u(t) T Mv(t) (14) 



(we are writing u,v for du/dt, dv/dt as is customary in mechanics); for a 
proof I refer you to your favorite calculus book. 

Let us now go back to the proof of the symplectic invariance of the first 
Poincare invariant. We write as usual the symplectic matrix S in block form 
A B\ 

q J so that the loop S~/(t) is parametrized by 

Wrt ( Ax{t) + Bp ^\ n<t<^ 
57(t) = [cx(t) + Dp(t) ) ' ° " * " ^ 

We thus have, by definition of the Poincare invariant, 

r-2n 

J(5 7 )= / (Cx(t) + Dp(t)f(Ax(t) + Bp(t))dt; 



Jo 

expanding the product in the integrand, we have I(S'y) = I\ + 12 where 

h = \ x(t) T C T Ax(t)dt + / p(t) T D T Bp(t)dt 
Jo Jo 

h = \ xitf C T Bp{t)dt + / p(tf D T Ax{t)dt. 
Jo Jo 



We claim that h = 0. Recall that C T A and C T B are symmetric in view of 
the two first equalities in (J7|); applying the differentiation formula (]14p with 
u = v = x we have 



[ 2n x(t) T C T A±(t)dt = \ ^{x{t) T C T Ax{t))dt 
Jo ^ Jo dt 



- [x(2ir)C T Ax(2ir) - x(0)C T Ax(0)} 



= 

because x(0) = x(2tt). Likewise, applying (fT4|l with u = v = p we get 

/•2-7T 
JO 



p(t)D T Bp(t)dt = 
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hence I\ = as claimed. We next consider the term I2. Rewriting the 
integrand of the second integral as 

x{t) T C T Bp(t) = p(t) T B T Cx(t) T 

(because it is a number, and hence equal to its own transpose!) we have 

-2tt r-2n 



I; 
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/ p{t) T B T Cx{tfdt + / p(t) T D T Ax(t)dt 
Jo Jo 



that is, since D T A = + B T C by transposition of the third equality in ([7]), 

i-2-k i-2tx 



2 







r Ztv 

p{t) T x{t)dt + / [p(t) T B T CAx(t) + p{t) T B T CAx{t)] dt. 
Jo 



Using again the rule (|14p and noting that the first integral is precisely 1(7) 
we get, D T A being symmetric, 



h = 1(7) + / * j t [ P (t) T B T CAx(t)] dt. 

The equality I(S'j) = 1(7) follows noting that the integral in the right-hand 
side is 

p{2tt) t B t CAx(27t) - p(0) T B T CAx(0) = 

since (x(2vr),p(2vr)) = (x(0),p(0)). 

The observant reader will have observed that we really needed all of 
the properties of a symplectic matrix contained in the set of conditions ([7]); 
this shows that the symplectic invariance of the first Poincare invariant is a 
characteristic property of symplectic matrices. 

2.3 Proof of the symplectic egg property 

Let us denote by Br the phase space ball centered at the origin and having 
radius R. It is the set of all points z = (x,p) such that \z\ 2 = |x| 2 + |p| 2 < R 2 . 
What we call a "symplectic egg" is the image S(Br) of Br by a symplectic 
matrix S. It is thus an ellipsoid in phase space, consisting of all points z 
such that S _1 z is in the ball Br, that is |S' _1 z| 2 < R 2 . Using formula (|10p 



giving the inverse of S = together with the relations A T C = C T A, 

B T D = D T B in d7J we get the following explicit expression after some easy 
calculations: 

x T (CC T + DD T )x - 2x T (DB T + CA T )p + p T (AA T + BB T )p < R 2 



9 



(don't worry: we will not have to use this cumbersome inequality in what 
follows!). 

Let us now cut S(Br) by a plane IX,- of conjugate coordinates Xj,pj. We 
get an elliptic surface Tj, whose boundary is an ellipse denoted by jj. Since 
that ellipse lies in the plane Hj we can parametrize it by only specifying 
coordinates Xj(t), Pj(t) all the other being identically zero; relabeling if 
necessary the coordinates we may as well assume that j = 1 so that the 
curve 7j can be parametrized as follows: 



for < t < 2vr with xi(0) = xi(2tt) and pi(0) = pi(2tt). Since x k (t) = 
and pk {t) = for k > 1 the area of the ellipse is given by the formula 



hence Area(Ti) = /(71). Since the inverse matrix S" 1 is symplectic, we 
have /(71) = I(S' _1 7i). But the loop S'~ 1 7i bounds a section of the ball Br 
by a plane (the plane S^Uj) passing through its center. This loop is thus 
a great circle of Br and the area of the surface S~ 1 Ti is thus exactly irR 2 , 
which was to be proven. 

We urge the reader to notice that the assumption that we are cutting 
S(Br) with a plane of conjugate coordinates is essential, because it is this 
assumption that allowed us to identify the area of the section with action. 
Here is a counterexample which shows that the property does not hold for 
arbitrary sections of S(Br). Take, for instance 



0; Pl (t),0,---,0) T 




Jo 

n p2-k 




71 
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1/Ai 
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, Ai > 0, A 2 > 0, and Ai / A 2 



(15) 



\0 



I/A2/ 



so that S(Br) is defined by the inequality 



^xj + — x\ + Aip? + \ 2 pl < R 2 . 



Ai A2 



10 



The section of S(Br) by the X2,P2 plane is the ellipse 

^x\ + Xxpl < R 2 
Ai 

which has area ir{R 2 yf\i -^/l/Ai) = -kR 2 as predicted, but its section with 
the X2,pi plane is the ellipse 

— x\ + \ 2 p\ < R 2 
Ai 

which has area ir(R 2 \J\\/\2) which is different from ttR 2 since Ai 7^ A 2 - 

The assumption that S is symplectic is also essential. Assume that we 
scramble the diagonal entries of the matrix S above in the following way: 



S' = 











° ^ 





A 2 














1/A a 





\o 











The matrix S' still has determinant one, but it is not symplectic (cf. (jlip ). 
The section S'(Br) by the X2,P2 plane is now the ellipse 

— x\ + \2P\ < R 2 
Ai 

with area ttR 2 ^ 'A1/A2 / ttR 2 . 

3 The Symplectic Camel 

The property of the symplectic camel is a generalization of the property of 
the symplectic egg for arbitrary canonical transformations; it reduces to the 
latter in the linear case. 

3.1 Gromov's non-squeezing theorem: static formulation 

As we mentioned in the Prologue, the property of the symplectic egg is 
related to the "non-squeezing theorem" of Gromov [T7] in 1985. To under- 
stand it fully we have to introduce the notion of canonical transformation 
[2, 4J. A canonical transformation is an invertible infinitely differentiable 
mapping 
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of phase space on itself whose inverse / 1 is also infinitely differentiable and 
such that its Jacobian matrix 

d(x',p') 



f'(x,p) 



d(x,p) 



is symplectic at every point (x,p). A symplectic matrix S = au- 
tomatically generates a linear canonical transformation by letting it act on 

phase space vectors ( X ) — > S I °° ) : it is an invertible transformation (be- 
W \PJ 

cause symplectic matrices are invertible), trivially infinitely differentiable, 
and the Jacobian matrix is the matrix S itself. Phase space translations, that 

is mappings ( °° ) — > ( x + x o\ j canonical: their Jacobian matrix is 
\pj \p + PoJ 

just the identity ^ ^ . By composing linear canonical transformations 

and translations one obtains the class of all affine canonical transformations. 

Here is an example of a nonlinear canonical transformation: assume that 
n = 1 and denote the phase space variables by r and (p instead of x and p; 
the transformation defined by (r,ip) — > (x,p) with 

x = V2r cos ip , p = V%r sin ip , < <p < 2tt, 

has Jacobian matrix 



f'(r, ip) 



= cos tp — h= sin ip 



2r sin ip \/2r cos ip 



which has determinant one for every choice of r and ip. The transformation 
/ is thus canonical, and can be extended without difficulty to the multi- 
dimensional case by associating a similar transformation to each pair (xj,pj). 
It is in fact a symplectic version of the usual passage to polar coordinates 
(the reader can verify that the latter is not canonical by calculating its 
Jacobian matrix); it can also be viewed as the simplest example of action- 
angle variable [21 H] ; for instance it reduces the isotropic harmonic oscillator 
Hamiltonian H = + x 2 ) to K = r. 

We will see in a moment why canonical transformations play such an 
important role in Physics (and especially in classical mechanics), but let us 
first state Gromov's theorem: 

Gromov's theorem: No canonical transformation can squeeze a ball Br 
through a circular hole in a plane lij of conjugate coordinates Xj,pj 
with smaller radius r < R . 
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This statement is surprisingly simple, and one can wonder why it took so 
long time to discover it. There are many possible answers. The most obvious 
is that all known proofs Gromov's theorem are extremely difficult, and make 
use of highly non-trivial techniques from various parts of pure mathematics, 
so the result cannot be easily derived from elementary principles. Another 
reason is that it seems, as we will discuss below, to contradict the common 
conception of Liouville's theorem, and was therefore unsuspected! 

So, what is the relation of Gromov's theorem with our symplectic eggs, 
and where does its nickname "principle of the symplectic camel" come from? 
The denomination apparently appeared for the first time in Arnol'd's paper 
[3]. Recalling that in [21] it is stated that 

'...Then Jesus said to his disciples, 'Amen, I say to you, it will be 
hard for one who is rich to enter the kingdom of heaven. Again 
I say to you, it is easier for a camel to pass through the eye of a 
needle than for one who is rich to enter the kingdom of God\ 

The biblical camel is here the ball Br, and the eye of the needle is the hole 
in the Xj,pj plane! (For alternative interpretations of the word "camel"; see 
the reader's comments following E. Samuel Reich's New Scientist paper [26] 
about [9].) 

Let us next show that the section property of the symplectic egg is indeed 
a linear (or affine) version of Gromov's theorem. It is equivalent to prove 
that no symplectic egg S(Br) with radius R larger than that, r, of the 
hole in the Xj,pj plane can be threaded through that hole. Passing S(Br) 
through the hole means that the section of the symplectic egg by the Xj,pj 
plane, which has area ttR 2 , is smaller than the area 7rr 2 of the hole; hence 
we must have R < r. 



3.2 Dynamical interpretation 

The reason for which canonical transformations play an essential role in 
Physics comes from the fact that Hamiltonian phase flows precisely consist 
of canonical transformations. Consider a particle with mass m moving along 
the x-axis under the action of a scalar potential V. The particle is subject 
to a force F = —-^;V(x). Since F = mdv/dt = dp/dt (Newton's second 
law), the equations of motion can be written 

dx dp dV 

dt dt dx 
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In terms of the Hamilton function 

H(x,p) = ^-p 2 + V(x) 

this system of differential equations is equivalent to Hamilton's equations of 
motion 

dx dH dp dH 

~di = ~dp~ ' dt = ~~dp~' ( ' 

We will more generally consider the n-dimensional version of (|17p which 
reads 

rim j fl TJ rlrt f) J-T 

1 < j < n. (18) 



dxj dH dpj dH 



dt dpj dt dxj 

(In mathematical treatments of Hamilton's equations [21 El IE] the function 
H can be of a very general type, and even depend on time t). In either case, 
these equations determine -as any system of differential equations does- a 
flow. By definition, the Hamiltonian flow is the infinite set of mappings <j)f 
defined as follows: suppose we solve the system (|18p after having chosen 
initial conditions xi(0), x n (0) and pi(0), ...,p n (0) at time t = for the 
position and momentum coordinates. Denote the initial vector thus defined 

x(oy 



, . Assuming that the solution to Hamilton's equations at time t 
p{0)J 

fx(t)\ 

exists (and is unique), we denote it by | |j j . By definition, qbf is just the 
mapping that takes the initial vector to the final vector: 



pit)) \p(0)J 



(19) 



As time varies, the initial point describes a curve in phase space; it is called 
a "flow curve" or a "Hamiltonian trajectory". 

The essential property to remember is that each mapping ct>^ is a canon- 
ical transformation; Hamiltonian flows are therefore volume preserving: this 
is Liouville's theorem [21 Uj. This easily follows from the fact that symplec- 
tic matrices have determinant one. Since it is not true that every matrix 
with determinant one is symplectic, as soon as n > 1 volume preservation 
also holds for other transformations, and is therefore not a characteristic 
property of Hamiltonian flows; see Arnold [2], Ch.3, §16 for a discussion of 
this fact. The thing to observe is that volume preservation does not imply 
conservation of shape, and one could therefore imagine that under the action 
of a Hamiltonian flow a subset of phase space can be stretched in all direc- 
tions, and eventually get very thinly spread out over huge regions of phase 
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space, so that the projections on any plane could a priori become arbitrary 
small after some (admittedly, very long) time t. In addition, one may very 
well envisage that the larger the number n of degrees of freedom, the more 
that spreading will occur since there are more directions in which the ball 
is likely to spread! This possibility, which is ruled out by the symplectic 
camel as we will explain below, has led to many philosophical speculations 
about Hamiltonian systems. For instance, in his 1989 book Roger Penrose 
([23]. p. 174-1 84) comes to the conclusion that phase space spreading sug- 
gests that " classical mechanics cannot actually be true of our world" (p. 183, 
1.-3). In fact, our discussion of Gromov's theorem shows that Hamiltonian 
evolution is much less disorderly than Penrose thought. Consider again our 
phase space ball Br. Its orthogonal projection (or "shadow") on any two- 
dimensional subspace II of phase space is a circular surface with area irR 2 . 
Suppose now that we move the ball Br using a Hamiltonian flow </>^ and 
choose for n the plane H, of conjugate coordinates Xj,pj. The ball will 
slowly get deformed, while keeping same volume. But, as a consequence of 
the principle of the symplectic camel, its "shadow" on any plane ILj will 
never decrease below its original value ttR 2 \ Why is it so? First, it is clear 
that if the area of the projection of /(Br) on on a plane Xj,Pj (/ a canonical 
transformation) will never be smaller than ttR 2 , then we cannot expect that 
/{Br) lies inside a cylinder (pj — aj) 2 + (xj — bj) 2 = r 2 if r < R. So is 
the "principle of the symplectic camel" stronger than Gromov's theorem? 
Not at all, it is equivalent to it! Let us prove this. We assume as in section 
12.31 that j = 1; this does not restrict the generality of the argument. Let 
71 be the boundary of the projection of /(Br) on the x\,p\ plane; it is a 
loop encircling a surface Y\ with area at least ttR 2 . That surface Tj can 
be deformed into a circle with same area using an area-preserving mapping 
of the xi,p% plane; call that mapping f\ and define a global phase space 
transformation / by the formula 

f(x 1 ,p 1 ,X 2 ,p 2 , ,X n ,p n ) = (fl(x 1 ,p 1 ),X 2 ,P2, ...,X n ,p n ). 

Calculating the Jacobian matrix it is easy to check that the matrix / is a 
canonical transformation, hence our claim. For a more detailed discussion 
of this and related topics see [9j [16] . 

4 Quantum Blobs 

By definition, a quantum blob is a symplectic egg with radius R = \ffi. 
The section of quantum blob by a plane of conjugate coordinates is thus 
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nh = \h. We will see that quantum blobs qualify as the smallest units 
of phase space allowed by the uncertainty principle of quantum mechanics. 
We begin with a very simple example illustrating the basic idea, which is 
that a closed (phase space) trajectory cannot be carried by an energy shell 
smaller (in a sense to be made precise) than a quantum blob. As simple as 
this example is, it allows us to recover the ground energy of the anisotropic 
quantum harmonic oscillator. 



4.1 The harmonic oscillator 

The fact that the ground energy level of a one-dimensional harmonic oscil- 
lator 

tt Px , 1 2 2 

H = — muj x 

2m 2 

is different from zero is heuristically justified in the physical literature by the 
following observation: since Heisenberg's uncertainty relation Ap^Ax > \% 
prevent us from assigning simultaneously a precise value to both position 
and momentum, the oscillator cannot be at rest. To show that the lowest 
energy has the value \hio predicted by quantum mechanics one can then 
argue as follows: since we cannot distinguish the origin (x = 0,p = 0) of 
phase space from a phase plane trajectory lying inside the double hyperbola 
p x x < , we must require that the points (x,p) of that trajectory are such 
that \p x x\ > \h\ multiplying both sides of the trivial inequality 

— + mux 2 > 2\px\ > h 
moj 



by w/2 we then get 



E = pL + -mu 2 x 2 > -hu 
2m 2 ~ 2 



which is the correct lower bound for the quantum energy. This argument 
can be reversed: since the lowest energy of an oscillator with frequency u 
and mass m is the minimal phase space trajectory will be the ellipse 



p 2 x + x 2 



mfojj {h/muj) 

which encloses a surface with area \h. Everything in this discussion imme- 
diately extends to the generalized anisotropic n-dimensional oscillator 

n P 2 1 

.7=1 3 
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and one concludes that the smallest possible trajectories in xj,pj space are 
the ellipses 

Pj , x l _ L 
rrijfkdj {h/rrijU]j) 

By the same argument as above, using each of the Heisenberg uncertainty 
relations 

ApjAxj > -h (20) 
we recover the correct ground energy level 

E = -FlWi + -hW2 H h -fiw n 

as predicted by standard quantum theory [22]. In addition, one finds that, 
the projection of the motion on any plane of conjugate variables Xj,pj will 
always enclose a surface having an area at least equal to |/t. In other 
words, the motions corresponding to the lowest possible energy must lie on 
a quantum blob! 



4.2 Quantum blobs and uncertainty 

The Heisenberg inequalities (|20p are a weak form of the quantum uncer- 
tainty principle; they are a particular case of the more accurate Robertson- 
Schrodinger [251 EZ] inequalities 

(Apj) 2 (Axj) 2 > A(xj,pj) 2 + jh 2 (21) 

(see Messiah [22] for a simple derivation). Here, in addition to the standard 
deviations Axj, Apj we have the covariances A(xj,Pj) which are a mea- 
surement of how much the two variables Xj,pj change together. (We take 
the opportunity to note that the interpretation of quantum uncertainty in 
terms of standard deviations goes back to Kennard [19]; Heisenberg's [IB] 
own interpretation was much more heuristic). Contrarily to what is often 
believed the Heisenberg inequalities (j20|) and the Robertson-Schrodinger in- 
equalities ([2T]) are not statements about the accuracy of our measurements; 
their derivation assumes on the contrary perfect instruments; see the dis- 
cussion in Peres |24j . p. 93. Their meaning is that if the same preparation 
procedure is repeated a large number of times on an ensemble of systems, 
and is followed by either by a measurement of Xj, or by a measurement of 
Pj, then the results obtained have standard deviations Axj, Apj; in addi- 
tion these measurements need not be independent: this is expressed by the 
statistical covariances A(xj,pj) appearing in the inequalities (|2TT) . 
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It turns out that quantum blobs can be used to give a purely geometric 
and intuitive idea of quantum uncertainty. Let us first consider the case 
n = 1, and define the covariance matrix by 

Ax 2 A(x,p) 
A(p, x) Ap 



2 )• (22) 



Its determinant is detS = (Ap) (Ax) — A(x,p) , so in this case the Robertson 
Schrodinger inequality is the same thing as det S > \h 2 . Now to the ge- 
ometric interpretation. In statistics it is customary to associate to £ the 
so-called covariance ellipse: it is the set of f^s points (x,p) in the phase 
plane satisfying 

hx,?)^ 1 ( X ) < 1. (23) 



2 V 1C ' \p 



Its area is 2tt\/ det S, that is 

Area(fi s ) = 2tt [(Ap) 2 (Ax) 2 - A(x,p) 2 ] 1/2 

and the inequality det S > \h 2 is thus equivalent to Area(i7s) > irh = \h. 
We have thus succeeded in expressing the rather complicated Robertson- 
Schrodinger inequality (I2ip in terms of the area of a certain ellipse. In higher 
dimensions the same argument applies, but contrarily to what common in- 
tuition suggests, the Robertson-Schrodinger inequalities are not expressed 
in terms of volume (which is the generalization of area to higher dimen- 
sions), but again in terms of areas -namely those of the intersections of the 
conjugate planes Xj,pj with the covariance ellipsoid 

'A(x,x) A(x,pY 



A(p,x) A(p, p) J ' 

Here A(x, x), A(x,p), etc. are the nx n block-matrices (A(xi, Xj)) lKi - <n , 
(A(xi,pj)) 1Ki j <n etc. Notice that the diagonal terms of E are just the 
variances Ax 2 ,...,Ax 2 l ;Ap 2 ,...,Ap 2 l so that (JM|) reduces to (i22l) for n = 1. 
Defining the covariance ellipsoid fix; as above, one then proves that the 
inequalities (|2~T|) are equivalent to the property that the intersection of fis 
with the planes at least \h. These inequalities are saturated (i.e. 

they become equalities) if and only if these intersections have exactly area 
T^h, that is, if and only if fis is a quantum blob! The proof goes as follows (for 
a detailed argument see [HI [16] ) : one first remarks, using a simple algebraic 
argument that the Robertson-Schrodinger inequalities are equivalent to the 
following condition of the covariance matrix, well-known in quantum optics, 
see e.g. [28, 29J and the references therein: 
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The eigenvalues of the Hermitian matrix £ + y J are non- 
negative: X + y J > 0. 

The next step consists in noting that in view of Sylvester's theorem from 
linear algebra that the leading principal minors of 

ih T _ ( A(x,x) A(x,p) + ^I 



S+ 2 J \A{p,x)-%I A(p,p) 

are non-negative. This applies in particular to the minors of order 2 so that 
we must have 



Ax] A(x J ,p J ) + § 

A(p v x 3 )-% Ap] 



> 



and this condition is precisely the Robertson-Schrodinger inequality (|2ip . 

As we have seen, the fact that the covariance ellipsoid is cut by the 
conjugate coordinate planes along ellipsoids with areas > \h implies the 
Robertson-Schrodinger inequalities. This is thus a geometric statement 
-and a strong one- of the quantum uncertainty principle, which can be 
rephrased as follows: 

Every quantum covariance ellipsoid contains a quantum blob, 
i.e. a symplectic egg with radius y/h. When this ellipsoid is a 
quantum blob, the Robertson-Schrodinger inequalities are satu- 
rated. 

This statement can be extended in various ways; in a very recent paper 
|13j we have applied this geometric approach to the quantum uncertainty 
principle to the study of partial saturation of the Robertson-Schrodinger 
inequalities for mixed quantum states. We show, in particular, that partial 
saturation corresponds to the case where some (but not all) planes of con- 
jugate coordinates cut the covariance ellipsoid along an ellipse with exactly 
area ^h; this allows us to characterize those states for which this occurs 
(they are generalized Gaussians). 

Another important thing we will unfortunately not be able to discuss 
in detail because of length limitations, is the following: everything we have 
said above still holds true if we replace the sentence "planes of conjugate co- 
ordinates Xj,pj" with the sentence "symplectic planes". A symplectic plane 
is a two-dimensional subspace of phase space which has the property that if 
we restrict the symplectic form to it, then we obtain a new symplectic form, 
defined on this two-dimensional space. For instance, it is easy to check that 
the symplectic planes (but those of coordinates xj,pk, j ^ k, or 



19 



Xj,X}~, Pj,Pk are not). One proves [21 |8] that every symplectic plane can 
be obtained from any of the Xj,Pj planes using a symplectic transforma- 
tion. This property implies, in particular, that the Robertson-Schrodinger 
inequalities (|2ip are covariant under symplectic transformations: if one de- 
fines new coordinates x',p' by (x',p') T = S(x,p) T , S a symplectic matrix, 
then if 

{A Pj ) 2 {Ax 3 ) 2 > A( X] , Pj ) 2 + lh 2 

we also have 

(a p ;.) 2 (ax;.) 2 > Aix'^p'^ + ih 2 . 

Also, there are possible non-trivial generalizations of the uncertainty 
principle, using new results in symplectic topology, for instance [T] which 
extends Gromov's theorem (in the linear case) to projections on symplectic 
subspaces with dimension greater than 2. In [12] we have shown how this 
result leads to "quantum universal invariants". 

5 Conclusion 

Quoting the great mathematician Hermann Weyl: 

l In these days the angel of topology and the devil of abstract 
algebra fight for the soul of each individual mathematical domain' 
(H. Weyl, 1939) 

This quotation from goes straight to the point, and applies to Physics 
as well: while algebra (in the large) has dominated the scene of quantum 
mechanics for a very long time (in fact, from its beginning), we are witnessing 
a slow but steady emergence of geometric ideas. Not only do these geometric 
ideas add clarity to many concepts, but they also lead to new insights (see 
e.g. |12j). This is what we had in mind while writing the present paper. 
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